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We derive an analytic solution for the ensemble-averaged collective dephasing dynamics of N noninteracting 
atoms in a fluctuating homogeneous external field. The obtained Kraus map is used to specify families of 
states whose entanglement properties are preserved at all times under arbitrary field orientations, even for states 
undergoing incoherent evolution. Our results apply to arbitrary spectral distributions of the field fluctuations. 
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Control of the coherent evolution of quantum systems in 
noisy environments [1] is one of the crucial prerequisites for 
exploiting nontrivial quantum effects in composite systems of 
increasing complexity. Whether in the context of controlled 
molecular reactions [2], of many-particle quantum dynamics 
[3], or of quantum computers and simulators [4], uncontrolled 
fluctuations and noise are detrimental to most purposes of op¬ 
timal control. Various strategies may be followed to coun¬ 
teract the harmful influence of the environment: shielding 
the system degrees of freedom [5], correcting environment- 
induced errors [6], exploiting basins of attraction in dissipa¬ 
tive systems [8], or compensating dissipation—e.g., by coher¬ 
ent dynamics [7], dynamical decoupling [9], or periodic mea¬ 
surements [10]. Such approaches can effectively reduce the 
environmental effects and can enhance coherence times, but a 
perfect protection of the quantities of interest is generally not 
possible. 

By restricting to superposition states within a decoherence- 
free subspace, initially entangled states can be shielded com¬ 
pletely from collective noise sources, hence protecting their 
entanglement at all times [11, 12]. However, such subspaces 
are rather fragile to small perturbations, which limit their ap¬ 
plicability in the context of dynamical processes [13]. Here, 
we identify conditions that ensure complete preservation of 
arbitrary degrees of entanglement, even for states that are not 
invariant under an incoherent time evolution. Specifically, we 
consider an important class of environment-induced fluctua¬ 
tions, which are frequently encountered in state-of-the-art ex¬ 
periments [14-16]: they manifest in intensity fluctuations of 
spatially homogeneous experimental control fields, giving rise 
to an effective dephasing process. We show how control of the 
external field’s orientation can lead to the complete preserva¬ 
tion of entanglement in bipartite—as well as multipartite— 
settings, for arbitrary spectral characteristics of the control 
field fluctuations. We further identify families of states ex¬ 
hibiting time-invariant entanglement for arbitrary orientations 
of the external field. 

To set the stage, let us consider a collection of N noninter¬ 
acting atomic two-level systems with identical energy split¬ 
ting Tiu> controlled, e.g., by a homogeneous magnetic field. 
Integration over the unavoidable fluctuations of the latter’s 
strength will induce a probability distribution p(to) of the 


characteristic energy splitting, and the /V-atom quantum state 
at time t therefore needs to be described by the statistical op¬ 
erator 

pit) = f p(oj)UJtf N p(0)Ul(tf N Aoj, ( 1 ) 

provided the field fluctuations occur on time scales which are 
longer than the time t over which the /V-atom state is prop¬ 
agated by the unitary f/ w (f) 0iv - In order to assess the open 
system time evolution of the quantum correlations inscribed 
into the /V-atom system, it is convenient to derive an explicit 
expression for p(t) in terms of the spectral distribution p(oj) 
characterizing the fluctuations. 

The single-atom propagator U,Jt) = e j s gener¬ 

ated by the time-independent single-atom Hamiltonian //„, = 
{fi(jj/2)n ■ cr, with cr = (cr x , cr y , cr.) being the vector of the 
Pauli matrices and n the orientation of the field. //„, de¬ 
scribes atomic dipoles interacting with electromagnetic fields, 
as, e.g., the electronic qubits in trapped-ion quantum reg¬ 
isters [14, 17], Introducing pairs of orthogonal projectors 
A ± = (I 2 ± n ■ cr)/ 2, we can rewrite the time evolution op¬ 
erator for a collection of N atoms as 

f7„(t) w = (e' i "' /2 A + + e ibJt,2 A_f N 

N 

s icv,U-N/2) & , > ( 2 ) 

j=0 

where we have defined the operators 

0 >= ( 3 ) 

J V J ’ S&.H 

where X, v denotes the symmetric group and V s = 

2,IT(i) ■ • • 4 (ao>< 4 • ■ ■ fjvl represents the permutation s e 
X)v in the operator space of N qubits. The ensemble-averaged 
state after time f, 

N 

pit) = J] M jk (t)&jp(0)& k , (4) 

j,k =0 

is then fully characterized by the Toeplitz matrix M(t), whose 
elements Mj k (t) = <p[{j - k)t] are generated by the character¬ 
istic function ip(t) — f p(co)e 1Mt dw of the probability distri¬ 
bution pica), Bochner’s theorem [18] ensures that M(t) is a 
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Hermitian semipositive definite matrix for all t. Diagonaliza- 
tion leads to the canonical Kraus form [19] 

N 

p(t) = 6,0 b(0)] = 2 A,(f)p(0)Aj( t ), (5) 

i=0 

where the Kraus operators A,-(f) = contain 

the eigenvalues A'(t) and the components of the eigenvectors 
A‘(t) of M(t). Note, from the structure of A,-(f), that the above 
defined Kraus operators mediate an effective interaction be¬ 
tween the individual qubits—with its origin in the spatial ho¬ 
mogeneity of the external field. These environment-induced 
interactions are able to create discord-type quantum correla¬ 
tions [17] and, as we will show in this Letter, given the appro¬ 
priate control of n, can uphold multipartite entanglement at all 
times for arbitrary intensity fluctuations. 

Using the fact that the operators A ± are orthogonal projec¬ 
tors in C 2 , we can immediately show that both of the oper¬ 
ators [<5)y}y and [A,(f)}, satisfy the condition 2,- Aj(t)A,-(t) = 
2/ ®j 0/ = Ij', which ensures that the map e,o defined in (5) — 
from now on called the “collective dephasing” map—is not 
only completely positive but also trace preserving for all t 
[19, 20], 

For absolutely continuous distribution functions [24], 
the characteristic function vanishes asymptotically; i.e., 
lim^co tp{t) = 0. We then have that lim,-,,*, M(t) = Ijv+i, and 
thus the Kraus operators reduce to lim^oo A,-(f) = 0,. The 
asymptotic A-qubit state is thus given by p s = lim^oopff) = 
Yji 0,p(O)©;. Because the operators 0, depend exclusively 
on the magnetic field direction n, the latter completely de¬ 
termines the properties of the asymptotic state. 

To gain some intuition on the time evolution of the entan¬ 
glement properties, e.g., of an A-ion quantum register un¬ 
der the action of the collective dephasing map, we first con¬ 
sider two-qubit states with maximally mixed reduced den¬ 
sity matrices (also called Bell-diagonal states). Such states 
allow for a simple geometric representation, since they are 
fully characterized by the matrix /3,y = tr(p • cr, ® crj) [1], 
There always exist unitary operations Ua and Ug such that 
UapU ] b has a diagonal /? matrix, /? = diag (di, c/ 2 , d$), while 
p and UapU b have the same separability properties [27], 
This allows us to associate with each density matrix a point 
d = (d ] ,d 2 ,d 3 ) r e R 3 and, because of positivity, any such 
point must lie inside a tetrahedron [Fig. 1 (a)] of vertices 
(-1,-1,-l) r , (-1,1, l) T , (1,-1, l) r and (1,1,-l) r , which 
represent the four Bell states [1], Inside this tetrahedron we 
distinguish an inner octahedron, which contains the separable 
states, from the four remaining corners, which consist of the 
entangled states [1] and are labeled by the Bell state they con¬ 
tain (e.g., T (-corner). In this setting Wootters’s concurrence 
[2] is simply the distance from the faces of the octahedron: 
C(d) = 1 /2 max(0, 2/ \di\ - 1 )■ Equidistant points, parallel 
to the surfaces of the octahedron, form the “isoconcurrence” 
planes. 

In the tetrahedron, the collective dephasing evolution is al¬ 
ways constrained onto a plane defined by Tr /3(t) = Tr/?(0) 


l® + > 



FIG. 1. Dynamics of a bipartite system with the initial state 
Po = I 4 / 2 O + 41®~) (®~| /5. (a) The evolution of the state un¬ 

der Lorentzian- (green) and box-distributed (red) noise, Co,i(w) and 
B 0 j(oj) (see the text for a definition), respectively, with a magnetic 
field direction n = (2, l,l) r /V6, is depicted inside the tetrahe¬ 
dron of Bell-diagonal states. The inner octahedron marks the set 
of separable states, (b) Depending on the orientation of the mag¬ 
netic field, the concurrence remains constant [for n = (1,0,0) r , or¬ 
ange], decays to a finite value [for n = (2, 1 , l) T / V6, blue], or de¬ 
cays to zero (for n = ( 0 , 1, l) T / V2, red). The decay is monotonic 
for Lorentzian (dashed lines), or nonmonotonic for box-distributed 
(solid lines) noise. 

[20], In the V I J (-corner, these planes coincide with isocon- 
currence planes, which implies that entanglement is preserved 
for all of these states, for arbitrary directions of the magnetic 
field. This leads to a finite-measure set of states with time- 
invariant concurrence, despite the fact that those states do 
evolve in time, p{t) 4 p(0) [29], For entangled states outside 
the |TL)-corner, we can use Eq. (4) to predict the final con¬ 
currence as Cd,f(n) = l/2max{0, 2 /Li0 ~ 2 n^)dj - 1], where 
n = 0 n\,n 2 ,n{) T and d = (d\,d 2 ,d^) T characterizes the ini¬ 
tial state [20], Thus, by solving for n, we can always find a 
field direction such that the entanglement is preserved at all 
times. This can be seen from the long-time limit in Fig. 1 (b), 
whereas the transient time evolution depends on p(«), as we 
will discuss later. 

Collective interactions become particularly relevant in mul¬ 
tipartite settings, where decoherence and dissipation can be 
strongly enhanced [11, 33-35], To analyze the effect of col¬ 
lective dephasing on multipartite entanglement, analytic ex¬ 
pressions a la Wootters [2] are not available. Intricate hier¬ 
archies of multipartite entanglement [36] can, however, be 
characterized efficiently by resorting to separability criteria 
















3 



FIG. 2. Influence of collective dephasing (5) on A-partite entangled 
IV-states. Here we show the upper bound k cff to the asymptotic state’s 
separability vs the polar angles 0 < 8 < tt/2 (measured from the z- 
axis to the x, y-plane) of the fluctuating magnetic field's direction. In 
blue and orange we show, respectively, the dependence of the critical 
angles 8 E and 0 npe, Eqs. (6) and (7), on A. The dots correspond 
to the numerical estimation of the smallest angle where Cfr changes 
from 1 to 2 (defining f? NPE ) or from IV - 1 to A (defining 8 E ). The 
uncertainty on this angle, due to the finite bin width of our sampling, 
is covered by the dot size. The lines represent the expressions (6) and 
(7). 


based on inequalities [37, 38], An A-partite state p is called k- 
separable if it can be written as a mixture of states of the form 
p = Pa, ®- ■ -®PA k , where A i ... Ak label a division of the A par¬ 
ties into k subgroups. For instance, the matrix elements in an 
arbitrary basis of any ^-separable A-qubit density matrix p sat¬ 
isfy 2 j0</< j<N—l Ip 2 '+l, 2 t+lI ^ Eo<i'</'<iV-l VFUP2 i +2-4l,2'+2-4l + 
(A - k)/2 E^o 1 P 2 - + i, 2 -+i [39, 40]. Defining k eff as the largest 
integer k saturating this inequality provides an upper bound 
to the states ^-separability class, as k < k e fj. When k e ff < 2 
the state certainly contains genuine multipartite entanglement, 
i.e., it is not even 2-separable, while the state can be fully sep¬ 
arable (A-separability) only if k e s > A. 

We consider the initial (A-partite entangled) VV'-state, IV) = 
([10.. .0) + |01.. .0) + ■ ■ ■ + |0.. .01))/ VA, where |1) and |0> 
denote eigenstates of <x.. Since the collective dephasing map 
(5) is invariant under the operation n —» -n and, additionally, 
this class of states exhibits rotational symmetry around the z- 
axis, the polar angle 6 e [0, /r/2] between n and the 2 -axis 
(which is defined by the local eigenbasis of the initial state) 
fully determines the evolution of the state under (5). Figure 2 
displays the entanglement properties of the resulting asymp¬ 
totic state, characterized by & e ff as a function of 6. In general, 
there are relatively small angle intervals that lead to a fully 
separable state, and typically k c a shows nonmonotonic depen¬ 
dence on 6. 

Our numerical data (Fig. 2) suggest that the asymptotic 
state resulting from | W) is certainly entangled (i.e. k clf < A) 
as long as 0 < 0e, where 

6e (A) = arctan (1 / Va) . (6) 


ArW k e ff(t) 



FIG. 3. (Color online) Evolution of the trace distance D tr between 
p(?) and the asymptotic state p s (solid lines), and of the state’s sepa¬ 
rability bound £ eff (connected dots), respectively, for eight-partite | W) 
(green) and | IV) (blue) initial states (see the text for a definition), box- 
distributed noise fluctuations Boj(ai), and a polar angle 8 = n/ 8 . The 
values where k eff > A, marked in a lighter green shade, indicate that 
the state is compatible with full separability, and hence they do not 
provide additional information than the case k eff = A. Note that the 
| IV) state exhibits time-invariant genuine multipartite entanglement 
outside of a time-invariant subspace. 


Conversely, when we choose a magnetic field that is close to 
the 2 -direction, the initial A-partite entanglement of the IV- 
state will be preserved during the dephasing process, since | W) 
is part of an eigenspace of the Hamiltonian for n = (0,0, l) r . 
Again, we find a critical angle 

#npe(A 0 = arctan (1 / y/N(N - 1)) , (7) 

such that for 8 < 0npe- the asymptotic state will contain gen¬ 
uine multipartite entanglement (£ e ff < 2). Conditions (6) and 
(7) provide a finite range of orientations that ensure preser¬ 
vation of entanglement properties in initial VV’-states. How¬ 
ever, as the number of qubits gets larger, higher accuracy is 
required to maintain A-partite entanglement (0 npe) or at least 
some type of entanglement (£?e)- Moreover, the fast decay 
of fl NPE with the number of qubits confirms that genuine A- 
partite entanglement is much more fragile than bipartite en¬ 
tanglement [36, 41], which is able to resist a larger range of 
field directions. We remark here that in order to modify 8 in a 
trapped-ion experiment it is much more natural to apply uni¬ 
tary pulses to the initial state to shift its relative orientation to 
the field, instead of actually changing the orientation of the 
external field [17]. 

Furthermore, we notice that states displaying time-invariant 
entanglement properties can be found in the multipartite case, 
too. One example is given by a specific family of IV- 
states, whose single-excited states carry the relative phases 
{e,' 2nk l N \^ =] in an arbitrary order, e.g., 

\W) =(e i(27r/JV) |10... 0) + e i(47r/A0 101... 0) 

+ + I0...01))/ VA. (8) 

As shown in Fig. 3, this state remains A-partite entangled 
throughout the whole evolution, but the state itself evolves 
into a stationary state, as is displayed by the trace distance 
Ar(0 = ||p(0 - PjII/ 2, where ||A|| = Tr VA^A denotes the 
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trace norm. The question remains whether this state is part 
of a finite-measure set of states whose multipartite entangle¬ 
ment properties are conserved, similarly to the l F )-corner 
in the bipartite case—notice that |W) reduces to 'F ) when 
N - 2. Such states would constitute ideal candidates for quan¬ 
tum computations by exhibiting invariance under collective 
dephasing effects. 

Let us finally characterize a family of time-invariant states, 
for arbitrarily many qubits. Using Eq. (4), it can be shown 
that any state of the form p w = £ veIv c s V s , where c s are arbi¬ 
trary coefficients and V s are the permutation operators defined 
above, satisfies p(t ) = p( 0) at all times [20]. These states, 
known as multipartite Werner states [3], are also characterized 
by their invariance under arbitrary local unitary transforma¬ 
tions U® N [42]. Since such transformations describe collec¬ 
tive changes of the local qubit coordinate systems, it is quite 
intuitive that these states are time-invariant for arbitrary direc¬ 
tions of the external field. This identifies a (N\ - l)-parameter 
family of states that always span a decoherence-free subspace 
[11, 12], In the geometric picture of Fig. 1 (a), these states lie 
on the line passing through the origin of the tetrahedron and 
the |TL)-state. 

We conclude with some remarks on the transient evolution 
towards the asymptotic state. To determine how close the 
evolved state is to its asymptotic state, we again employ the 
trace distance Z) tr (f) which has a clear interpretation in terms 
of the distinguishability of the quantum states [43], In our 
present context, the trace distance is employed as an autocor¬ 
relation function, which reveals the monotonicity of the quan¬ 
tum evolution. 

While different types of noise fluctuations lead to the same 
asymptotic state, as discussed earlier, the transient behavior 
can be qualitatively different, as displayed in Fig. 1. When the 
distribution p(u >) is Lorentzian, C WOir (w) = (y/n)[(co - a> o) 2 + 
y 2 ] -1 , or Gaussian, N t Jo>cr (Y) = exp(-(w - a>o) 2 /2cr 1 )/ V2;rcr 2 
(as suggested in Ref. [41]), the properties of the state, such 
as the concurrence, decay exponentially towards their asymp¬ 
totic value [Fig. 1 (b)]. When we instead consider the 
box distribution over the interval [0, o>o], he., Bo,w 0 ( w ) = 
[@//(w) - ©h(w - wo)] /wo, where ©h(w) is the Heaviside 
step function, we observe a nonmonotonic approach of the 
quantum system to the asymptotic state (Fig. 1). In fact, for 
this distribution the characteristic function i p(t) is proportional 
to sin(woO/L which asymptotically decreases on a signifi¬ 
cantly longer time scale than the exponential decay charac¬ 
terizing the Lorentzian or Gaussian distributions. This non¬ 
monotonic behavior also implies that the ensemble-averaged 
dynamics of noninteracting atoms in a fluctuating classical 
field cannot be modeled by an effective Markovian environ¬ 
ment for certain noise distributions p(co) [34, 44]. These fre¬ 
quency fluctuations, therefore, take on the role of the environ¬ 
ment’s spectral density in a standard open-system description 
of decoherence [34, 35], 

To summarize, we have provided a model for the dephasing 
dynamics of a collection of noninteracting atoms subject to a 
homogeneous external field of fluctuating intensity. The effec¬ 


tive environment-induced interactions are described analyti¬ 
cally by an exact solution in terms of a canonical Kraus map, 
able to describe the time evolution of multipartite systems un¬ 
der arbitrary intensity fluctuations. Our model applies to a 
variety of experiments in atomic physics, and describes one 
of the dominant error sources for state-of-the-art trapped-ion 
experiments. Complete theoretical control on transient as well 
as asymptotic dynamics allows for the formulation of precise 
conditions for preserving relevant quantities, such as entan¬ 
glement, as well as for the identification of families of states 
whose properties are completely insensitive to the direction of 
the external field. 
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Supplementary Material 

TRACE PRESERVING PROPERTY OF THE COLLECTIVE 
DEPHASING MAP 


To prove that tr P(t) is conserved under time evolution, we 
compute its derivative and check whether it vanishes: 


dtr/Xt) 

df 


~4 tr [p(t) |'T / ) < V F |]. 


(14) 


The operators A ± = j(I 2 ± n • tr) form a complete set of 
orthogonal projectors on the Hilbert space 'H - C 2 of each 
qubit. These properties are inherited by the 0, operators, 
which are themselves orthogonal projectors: 0 , = 0 1 and 

0/0 j - 

The trace preserving property of the map therefore reduces 
to 2 ; 0 ; 0 ; = 21 0 i = This is simply proven by putting 
t = 0 in Eq. (2) from the main text, and following the equali¬ 
ties from left to right: 

N 

2©; = (A + +A_) 0A, = I 2 *. (9) 

(=0 

For the A,-(f) operators we instead have 

N _ 

YjKiOAM = X X 4 !( f)T‘ ( t)[4 ( r)]*0,0 ; 

1=0 i j,k 

= X X 

i j 

= ^0j = V, (10) 


The time derivative of p(t) reads 

N 

pit) = z, \pit)} = ^ Mijimpmj , d5) 

U =0 


which implies that 


= tr 


= tr 


tr[ J 6 ( 0 |T'-)<'n] 

N 

^%) 0 ,p(O 0 ,|T-)(r| 
ij =0 

N 

Y M.jiDQj |T"> (T~| 0,p(t) 

,j=o 

= tr[( 2 f [|T'-)<T'-|]) t p( 0 ]. 


(16) 


However, the Bell state T ) is an eigenstate of the Hamilto¬ 
nian of the system for whichever choice of n, which means 
that it is itself unaffected by collective dephasing. This im¬ 
plies that Zt | IT ) (T |] = 0, Vr and therefore 


dr 


tr p(t) = 0 . 


(17) 


where we have used the spectral decomposition M(t) = 
2; Z(t)A'(t)A^(t), together with Mjj(t) = 1 ,Vf. 


CONSERVED TRACE OF THE /? MATRIX 


We now prove that the trace of the /3it) matrix, defined by 
Pi jit) = tr [pit) ■ <Tj <S> (T j |, is a time-invariant quantity. From 
the definition we have 


tr Pit) = ^]/?;;(0 = Y tr W0 • (T i ® <rd 

i= 1 1=1 

3 

pit) Y^ CTi ® <- 


= tr 


i= 1 


(ii) 


Notice now that the Bell state IT ) (T | reads [1] 

3 


|r-><r-| = k 


-2' 


CT; S) CTj 
(=1 ) 


( 12 ) 


which yields 2; (T i ® cr,- = I 4 — 4 |T ) (T |. Substituting back 
we then have 


ti Pit) = tr pit)-4tr [pit) ]*¥-) (*¥-[] 

= 1 -4trLp(0|T ><T-|]. (13) 


CONCURRENCE FOR BELL-DIAGONAL STATES 


In the tetrahedron of Bell-diagonal states [1], Wotters’ con¬ 
currence [ 2 ] reads 


dd)= ^ max jo, -l+^k/elf, 


(18) 


where d = id\,didZ) T e R . 3 is the point representing the 
quantum state (see main text). 

For states in the T (-corner, we have d, < 0, V;. Hence, the 
concurrence in this corner can be rewritten as 

Cid) = ~ max jo, -1 - ^ dj 

= -max {0,-1 -k), (19) 

where we have used the fact that tr /? = 2 ; dj — k is a constant. 
This explicitly proves that Bell-diagonal states in this corner 
have time-invariant concurrence. 

In the other corners of the tetrahedron, only one of the co¬ 
ordinates is negative. If we suppose that d 1 < 0 (i.e. |<f>_)- 
corner), we have 

— 1 + ^ \di\ — —1 — d\ +r/ 2 +t/3 = —1 + k — 2 d\. 
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Let d represent the initial state, whereas d' represents the final, 
asymptotic state of the system after collective dephasing. Fur¬ 
thermore, we denote the negative components of d and d with 
a subscript j, i.e., we have dj < 0 and A < 0, respectively. Di¬ 
rect application of the collective dephasing map leads to 



( 20 ) 


The concurrence in the final state is then 


C (</ f ) = max 10, -1 + y dj - 2 y d^ 


- ^maxjo,-l + ^(1 -2nf)dj\. ( 21 ) 


From these formulae one can immediately notice that C(d) = 
C(d { ) <=> n = ±e j, where {^1,62,63) is the standard basis 


of R 3 . 


TIME-INVARIANCE OF MULTIPARTITE WERNER 
STATES 


In this Section we prove that the multipartite Werner states 
[3] are time-invariant under the action of the collective de¬ 
phasing map. Let s e Zjv be a permutation in the symmetric 
group, and V s its representation in the operator space of N 
qubits: 


V s - ^ |i s (i) • • ■ i s (.N)Xh ■ ■ ■ Ovl- (22) 

i| ...!A-6{0,1} 

The multipartite Werner states are then defined as pw — 
Y A ci v c s V s , where c s are arbitrary coefficients, leading to a 
valid quantum state p w . 


Analogously, if we define k, = i\(N - i)l and <3; = Af' ® 
we can rewrite the 0, operators as 

0 / - l X V ' Q ‘ V >- (23) 

* S&N 

Direct application of the collective dephasing map yields 

N 

e[p w ] = Y, Z ^{ v ^ v l) v 4 v ^ Q ‘ v D- (24) 

i=0 n,<T,AeH^ / 

Because X.y is a closed group, the concatenation of two per¬ 
mutations describes another permutation, and therefore 3a e 
■ Vcr — V„V„. The expression above can then be rewritten 
as 


e\Pw\ 


= £ 2 %M v A v D( v & v 'k) 

1=0 7r,a,Aeh^ i 

‘Y c - V -YY 

nel ,n i—0 cheLm 

N 

~ ^ ^ n ^ ^ — Pw 9 


(K&vI) (VaQV*) 


h • ' ^ Ic 

i=0 ael,v ' /teljv 1 


(25) 




i=0 


where we have used the idempotency of the 0, operators and 
the closure relation Yi = L ;V • 
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